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ABSTRACT 

We use the outer gap model to explain the spectrum and the energy dependent light curves 
of the X-ray and soft 7-ray radiations of the spin-down powered pulsar PSR B1509-58. In the 
outer gap model, most pairs inside the gap are created around the null charge surface and the 
gap's electric field separates the two charges to move in opposite directions. Consequently, the 
region from the null charge surface to the light cylinder is dominated by the outflow current while 
that from the null charge surface to the star is dominated by the inflow current. We suggest that 
the viewing angle of PSR B1509-58 only receives the inflow radiation. The incoming curvature 
photons are converted to pairs by the strong magnetic field of the star. The X-rays and soft 7-rays 
of PSR B1509-58 result from the synchrotron radiation of these pairs. Magnetic pair creation 
requires a large pitch angle, which makes the pulse profile of the synchrotron radiation distinct 
from that of the curvature radiation. We carefully trace the pulse profiles of the synchrotron 
radiation with different pitch angles. We find that the differences between the light curves of 
different energy bands are due to the different pitch angles of the secondary pairs, and that the 
second peak appearing aX E > lOMeV comes from the region near the star, where the stronger 
magnetic field allows pair creation to happen with a smaller pitch angle. 

Subject tieadings: methods: numerical -pulsars:general -radiation mechanisms; non-thermal 



1. Introduction 

PSR B1509-58 (hereafter PSR B1509) is a 
unique subject of the multi-wavelength study of 
the high energy radiation mechanism of the spin- 
down powered pulsar. It was discovered in soft 
X-ray band by Einstein (Seward & Harnden 1982) 
and soon was detected in the radio band (Manch- 
ester et al. 1982) thirty years ago. PSR B1509 is 
a young (about 1600 years old) spin-down pow- 
ered pulsar with a period of about 150 ms and a 
surface magnetic field of ~ 1.5 x lO^^G, which is 
much stronger than those of most of the canoni- 
cal pulsars (^ lO^^G). After being discovered, its 
radiation was observed by many X-ray and soft 
7-ray telescopes: EXOSAT (Trussoni et al. 1990), 
Ginga (Kawai et al. 1992), BASTE and OSSE 
on the CGRO (Ulmer et al. 1993; Matz et al. 
1994), SIGMA (Laurent et al. 1994), Welcome-1 



(Gunji et al. 1994), ROSAT (Greiveldinger et al. 
1995), RXTE (Marsden et al. 1997; Rots et al. 
1998; Livingstone & Kaspi 2011), ASCA (Saito et 
al. 1997), COMPTEL (Kuiper et al. 1999) and 
BeppoSAX (Cusumano et al. 2001). These obser- 
vations show that, from 0.1 keV to 10 MeV, PSR 
B1509 has a wide single peak, whose shape can be 
described by two Gaussian functions, and COMP- 
TEL shows that, when E > lOMeV, another peak 
shows up near the peak of radio emission. Af- 
ter 2008, the bottleneck in 7-ray astrophysics was 
broken by Fermi and AGILE. In 2010, these two 
7-ray telescopes provided the observed data of 
PSR B1509 above 30MeV (Abdo et al. 2010; Piha 
et al. 2010). The soft 7-ray radiation and the 
second peak found by COMPTEL have been con- 
firmed by them. 

PSR B1509 is unique among the known 7- 
ray spin-down powered pulsars for its very soft 
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7-ray spectrum and a very strong surface mag- 
netic field. It is suggested tliat tliere may be a 
new class of spin-down powered pulsar with high 
magnetic field (>10"'^^G) and soft 7-ray radiation, 
called soft 7-ray pulsar (Pilia & Pellizzoni 2011). 
Also, there is at least another pulsar that is simi- 
lar to PSR B1509, PSR J1846-0258, a young pul- 
sar with P ~ 324ms and surface magnetic field 
Bg ~ 4.9 X lO^^G, which cannot be detected by 
Fermi but is seen by X-ray telescopes (Kuipcr & 
Hermsen 2009; Parent et al. 2011). Therefore, 
a multi-wavelength study of the radiation mecha- 
nism of PSR B1509 can be a bridgehead to inves- 
tigate this possible new class of pulsars. 

Recent observations show that two kinds of pul- 
sars, the spin-down powered pulsar with Bg ~ 
lO^^G and the magnetar with B^ ~ 10" - lO^^G, 
have no clear boundary between them. It is be- 
lieved that the radiation of a magnetar (AXP or 
SGR) is powered by the energy stored in the su- 
per strong magnetic field with Bg > lO^^G (Dun- 
can & Thompson 1992), however, two low mag- 
netic field soft 7-ray repeaters, SGR 04184-5729 
with Bs < 7.5 X lO^^G (Rea et al. 2010) and 
J1822. 3-1606 with B,, ~ 2.7 x lO^'^G (Rea et al. 
2012), have been found, which challenge the re- 
quirement of a super strong magnetic field of the 
existing theoretical models. Furthermore, radio 
emission, which is usually regarded as a charac- 
teristic behavior of a spin-down powered pulsar, 
has been found in three magnetars, XTE J1810- 
197 (Camilo et al. 2006), IE 1547-5408 (Camilo et 
al. 2007) and PSR J1622-4950 (Levin et al. 2010). 
On the other hand, PSR J1846-0258, a spin-down 
powered pulsar with Bg ~ 4.9 x lO^^^G, is found to 
have a magnetar-like outburst (Kuiper & Hermsen 
2009). So, the study of the radiation mechanism 
of the very high magnetic field spin-down powered 
pulsars, and the connection between the canoni- 
cal spin-down powered pulsar and the magnetar, 
can help us to understand the physics of the two 
groups of neutron star. As a family member of 
the high magnetic field spin-down powered pul- 
sars, PSR B1509 may provide some insights for 
understanding the differences between typical pul- 
sars and magnetars. 

As a known 7-ray pulsar before the Fermi era, 
PSR B1509 has been interpreted by two main ra- 
diation scenarios, the polar cap model (Ruderman 
k Sutherland 1975; Daugherty & Harding 1982) 



and the outer gap model (Cheng, Ho & Ruderman 
1986; Romani 1996; Hirotani 2006). For the po- 
lar cap model, where the 7-ray photons are emit- 
ted near the stellar surface, the soft 7-ray pho- 
tons are generated via the photon splitting process 
7 — >■ 77 caused by a strong magnetic field (Hard- 
ing & Baring 1997). In the point of view of the 
outer gap model, where the 7-ray radiation comes 
from the outer magnetosphere, Zhang and Cheng 
(2000) assumed PSR B1509 to be another Crab 
pulsar and explained the X-ray to soft 7-ray radi- 
ation of this pulsar as the synchrotron radiation of 
the secondary pairs created by photon-photon pair 
creation. Since the light cylinder of PSR B1509 is 
much larger than that of the Crab pulsar, the in- 
verse Compton component of PSR B1509 is much 
weaker than that of the Crab. They predicted that 
the GeV energy flux produced by inverse Comp- 
ton scattering is nearly three orders of magnitude 
weaker than that of the MeV energy flux. However 
this predicted energy flux is nearly two orders of 
magnitude lower than the recent observations by 
AGILE (Pilia et al. 2010) and Fermi (Abdo et al. 
2010). It is interesting to note that the observed 
total radiation power of PSR B1509 is less than 
1% of its spin-down power, which is much lower 
than the efflciency predicted by most theoretical 
models (e.g. Zhang & Cheng 1997; Takata, Wang 
& Cheng 2010). Furthermore, the light cylinder 
of PSR B1509 is so large as to make it very diffi- 
cult to attenuate all of the GeV curvature photons 
emitted from the outer gap. 

Cheng, Ruderman & Zhang (2000) showed that 
most electron/positron pairs created inside the 
outer gap are produced near the null charge sur- 
face, the strong electric field inside the gap sepa- 
rates the two charges to move in opposite direc- 
tions. Therefore from the null charge surface to 
the light cylinder the radiation is mainly outward, 
and from the null charge surface to the star the ra- 
diation is mainly inward. Since the light cylinder 
of PSR B1509 is much bigger than that of the Crab 
pulsar, most outflow curvature photons can avoid 
photon-photon pair creation. If the line of sight 
is in the direction of the outgoing radiation beam, 
the spectrum of PSR B1509 should be a charac- 
teristic pulsar spectrum, namely a power law with 
exponential cut-off as predicted by Wang, Takata 
& Cheng (2010). Obviously this conflicts with the 
observed properties of PSR B1509. In this paper 
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we propose that the viewing angle of PSR B1509 
is in the direction of the incoming radiation beam, 
where the curvature photons with hundreds MeV 
are converted into pairs by the strong magnetic 
field. The observed X-rays and soft 7-rays are the 
synchrotron radiation of these pairs. It is also im- 
portant to note that the electric field from the null 
charge surface decreases rapidly towards the in- 
ner boundary (e.g. Hirotani et al. 2003; Takata, 
Chang & Cheng 2007; Tang et al. 2008), therefore 
the incoming radiation flux is at least an order of 
magnitude weaker than that of the outgoing flux. 

Most of these incoming curvature photons emit- 
ted by canonical pulsars with a surface magnetic 
field ~ lO^^G will be converted into pairs and 
hence they become invisible. On the other hand 
if they are not absorbed they may appear near the 
phase position of the radio pulse. Figure |3] gives 
the pulse phases and viewing angles of the incom- 
ing curvature photons, which are represented by 
the gray lines. The pulse phase of the center of the 
polar cap is set to be 0, because the radio emis- 
sion is believed to come from the polar cap. If the 
viewing angle is larger than the inclination angle of 
the pulsar and there is no magnetic pair creation, 
the gray lines will make a peak near the phase 0. 
There are some clear evidences for the existence 
of this incoming curvature photon beam in the en- 
ergy range of a few hundred MeV in gamma-ray 
millisecond pulsars whose surface magnetic field is 
of order of Bs ~ 10**G (Abdo et al. 2010b) but 
this mangetic field strength is not strong enough 
to convert photons of several hundred MeV into 
pairs. Even very strong multiple fields may exist 
on the surface of millisecond pulsars but they de- 
crease much more rapidly than the dipolar field. 
The incoming photons of several hundred MeV can 
escape from the magnetic pair creation process un- 
less they come closer than a distance of 2-3 stellar 
radu (Takata, Wang & Cheng 2010). 

If the pairs are generated by photon-photon 
pair creation, which is used to explain the radi- 
ation of the Crab pulsar in Tang et al. (2008), the 
pulse phases of the peaks of X-ray radiation are 
the same as those of the 7-ray radiation, because 
of the small pitch angles of the secondary pairs. 
On the other hand, the magnetic pair creation re- 
quires the photon to penetrate with a large angle 
to the magnetic field, therefore the incoming 7- 
ray photons have to go inwards to obtain enough 




Fig. 1. — The sketch of the structure of the gap. 
Outside the null charge surface the structure of 
the gap can be simplified as a two-layer structure. 
From the null charge surface to the stellar surface, 
the accelerating electric field is screened out. 

B±, and the secondary pairs obtain larger pitch 
angles than those from photon-photon pair cre- 
ation. These pitch angles make the pulse profile of 
the synchrotron radiation of these secondary pairs 
different from that of the curvature radiation of 
the primary particles. It has been found that the 
non-thermal X-ray pulse profile of the spin-down 
powered pulsar can be different from that of the 
7-ray, such as PSR J0659-I-1414 and PSR J1420- 
6048 (Weltevrede et al. 2010). This indicates that 
the radiation mechanism discussed in this paper 
can be further applied to explain the non-thermal 
X-ray radiation of other 7-ray pulsars. 

In section 2, we describe our theoretical model. 
Section 3 gives the details of the calculation of the 
pulse profile of the synchrotron radiation emitted 
by the secondary pairs generated by the magnetic 
pair creation process. In section 4, we present our 
calculation of the spectrum and the energy depen- 
dent light curves of PSR B1509. We discuss some 
special features in the observed spectrum and light 
curves in section 5. A brief summary is given in 
section 6. 

2. Theoretical Model 

2.1. The structure of the gap 

Figure [T] illustrates the structure of our the- 
oretical model. From the stellar surface to the 
light cylinder, along the magnetic field line, the 
emission region has different structures on the two 
sides of the null charge surface. For PSR B1509, 
only the radiation from the region inside the null 
charge surface is received by the viewing angle 
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while that from the region outside the nuh charge 
surface is missed. Outside the null charge sur- 
face, the particles go outwards along the open field 
lines, and the structure of the gap can be simpli- 
fied as a two-layer structure (Wang et al. 2010, 
2011) which contains a main acceleration region 
and a screening region. The main acceleration re- 
gion, where the charge density is much lower than 
the Goldreich- Julian charge density, pcj ^ ~Tnc 
(Goldreich & Julian 1969), is the source of the 
curvature photons with E > IGeV, and a strong 
accelerating electric field is formed from such de- 
viation. Above the main acceleration region, the 
screening region has a charge density much larger 
than pgj, which screens out the electric field in the 
trans-field direction. In this region, the charged 
particles near the edge between the two layers can 
be accelerated to emit curvature photons of hun- 
dreds of MeV. Inside the null charge surface, the 
electric field is reduced and there is no two-layer 
structure. 

In the direction along the magnetic field line, 
the accelerating electric field is screened out in a 
similar way as in the trans-field direction outside 
the null charge surface. The potential that gen- 
erates the accelerating electric field in the gap is 
given by the Poisson equation 'V^^' = — 47r(p — 
Pgj)- In the region of two-layer structure, ac- 
cording to this Poisson equation, in the trans- 
field direction, a strong accelerating electric field 
is formed in the region where the charge density is 
close to zero (the main acceleration region), and 
screened out in the region where the charge den- 
sity is larger than poj (the screening region). The 
average charge density of the two-layer region is 
~ 0.5pGj (Wang et al. 2010). In the same way, 
in the direction along the field line, the accelerat- 
ing electric field is formed outside the null charge 
surface, where p ~ O.bpGj, and screened out in- 
side the null charge surface where |p| increases. It 
is still unclear what the distribution of the charge 
density and the electric field below the null charge 
surface are, and how the region with the two-layer 
structure connects to the region below the null 
charge surface. However, in this paper we discuss 
the radiation processes far from the null charge 
surface, so we can leave such questions to a future 
general study of the soft 7-ray and X-ray radia- 
tions from the spin-down powered pulsars. 

The origin of these charged particles between 



the stellar surface and the null charge surface has 
been discussed by Takata et al. (2010). The 
photon-photon pair creation between the thermal 
X-ray photons and the GeV photons generates 
pairs near the null charge surface, and the mag- 
netic field generates pairs near the stellar surface 
via magnetic pair creation. The magnetic field 
of the pulsar has two parts: the dipole field and 
the multi-pole field. The dipole field, which can 
be determined from the observed period P and 
its derivative P, dominates the magnetosphere of 
the neutron star, and determines the accelerating 
electric field together with the period of the pulsar 
P. The multi-pole field, which cannot be observed 
directly, dominates the region near the stellar sur- 
face and plays an important role in the gap closure. 
The charged particles under the null charge sur- 
face move inward to the stellar surface along the 
magnetic field lines, losing their energy by the cur- 
vature radiation. Because the accelerating electric 
field is weak, the lost energy cannot be compen- 
sated for by the electric field, and the energy of 
the incoming curvature photons emitted far from 
the null charge surface is around hundreds of MeV. 
If these curvature photons are converted to pairs 
by the multi-pole field, part of the secondary pairs 
can obtain outgoing momentums from the multi- 
pole field and go outwards along the field lines 
to close the gap. These pairs may radiate syn- 
chrotron photons after being generated, but they 
can be neglected in this model, because 1) the pho- 
tons may be converted to pairs by the multi-pole 
field, 2) they may be covered by the neutron star, 
and 3) the pairs spiral around the complex multi- 
pole field lines. These synchrotron photons cannot 
show the peak caused by the caustic effect, which 
is due to the photons emitted from different places 
in the magnetosphere having parallel directions. 

However, not all of the incoming curvature pho- 
tons arc converted to pairs by the multi-pole field; 
if the dipole field is strong enough the curvature 
photons can also be attenuated by the dipole field. 
If the pairs arc generated by the dipole field, their 
synchrotron radiation has a chance to avoid enter- 
ing the region dominated by the multi-pole field, 
and shows a significant pulse profile due to the 
dipole field's caustic effect. 
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2.2. The magnetic field and the electric 
field 

We adopt the three dimensional rotating vac- 
uum dipole field as the global magnetic field, 
which is given hy B ^ r[f ■ + §^ + ■^)] - 
(^ + ^ + J?) (Cheng, Rudermanfe ZhangVoOO), 
where fj, = ij,{x sin a cos flt + y sin a sin flt + z cos a) 
is the magnetic moment vector, f is the radial unit 
vector, and a is the inclination angle. We do not 
know the relationship between the surface multi- 
pole field and the global dipole field. For PSR 
B1509, a pulsar with a strong dipole field, we as- 
sume the multi-pole field dominates the region of 
r < 1.5i?s, where r is the distance to the center of 
the neutron star, and Rg is the stellar radius. 

The realistic magnetic field of the pulsar's mag- 
netosphere is a compromise between the force 
free field and the rotating vacuum dipole field, 
which are two eigenstates of the magnetic field 
of the spin-down powered pulsar. Numerical re- 
sults show that the two fields have distinct field 
structures, which provide different pulse profiles 
(Spitkovsky 2006; Bai & Spitkovsky 2010). How- 
ever, our model of the synchrotron radiation does 
not strongly require a realistic field structure, be- 
cause the magnetic pair creation and the syn- 
chrotron radiation happen near the star, where the 
structures of the force free field and the vacuum 
dipole are close to each other. 

For the accelerating electric field outside the 
null charge surface, we apply our three dimen- 
sional two-layer structure model (Wang et al. 
2011). Although the outgoing curvature photons 
of PSR B1509 are missed by the viewing angle, 
for the integrity of the model we here repeat the 
equations: 
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for hi < z < h2 
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The X, z and (j)p represent coordinates along the 
magnetic field line, the height measured from the 



last-open field line, and the azimuthal direction. 
hi and /i2 represent the thickness of the main- 
accretion region and the total gap thickness, re- 
spectively. The g = [p— Pgj) I Pgj is the deviation 
of the charge density from the Goldreich- Julian 
charge density pcj and in the unit of pcj'- 



g(z,4>p) = 



-51 (0p) 
.92 (0p), 



if 
if 



{)<z<hi 
hi < z < h2 



(2) 



where gi > 0, g2 > and they satisfy (ft.2/^1)^ — 
1 + 51 /(?2- The /i2 and the polar radius rp satisfy 
the relationship, / = h2{Rs)/rp. 

We do not consider the distributions of the 
three parameters, 1 — gi, hi/h2 and /, in the az- 
imuthal direction in the gap, which are found from 
the energy dependent light curves of the Vela pul- 
sar (Wang et al. 2011). This is because the soft 
7-ray photons of PSR B1509 are not the outgo- 
ing curvature photons. We use 1 — gi — 0.05, 
hi/h2 = 0.92, and / = 0.35. 

Under the null charge surface, there is so far 
no picture about the distribution of the charge 
density and the corresponding reduced accelerat- 
ing electric field as clear as the two-layer structure 
for the region above the null charge surface. The 
electric field is assumed to decrease quadratically 
along the field line to the inner boundary of the 
gap, where i^y = 0. We adopt the form of the elec- 
tric field used in Tang et al. (2008) and Takata et 
al. (2010): 



E\\ (r < r„ 



1 



{rnull/rinY - 1 



E\\^nulU (3) 



where the and r„„n are the distance of the in- 
ner boundary of the gap and the null charge sur- 
face from the star, respectively, and the is 
the electric field at the null charge surface, which 
is given by Equation ([T|). For PSR B1509, which 
has a period of 150ms, the position of the inner 
boundary is not important here, because the en- 
ergy of the curvature photons emitted by the in- 
coming particles are not determined by the weak 
electric field in the large magnetosphere, except in 
the region close to the null charge surface. We set 
the inner boundary as = 20i?s. 

The electric field defined above can still accel- 
erate the particles near the null charge surface to 
emit curvature photons even up to >lGeV. These 
particles go towards the stellar surface along the 
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field lines and lose their energy through emitting 
curvature photons, whose energy drop down from 
GeV to hundreds of MeV. In order to know the 
number of these particles at any place between 
the stellar surface and the null charge surface, we 
set the number density of the incoming particles 
as, 

n(r<r„„„)=77— (4) 
Znec 

which is of the same order as the local Goldreich- 
Julian charge density, and i] is a fitting parame- 
ter. Because of the conservation of the magnetic 
flux, the flux of these incoming particles is a con- 
stant along the field line, from which we can cal- 
culate their curvature radiation. Along one field 
line, 77 should be continuous at the null charge sur- 
face, because the outgoing and incoming particles 
are generated from the pair creations occurring 
around the null charge surface. If we define the 
charge density above the null charge surface as 
positive, it nearly equals the density of the out- 
going positrons, which can be described as p{r > 
Tnuii) ~ P+- On the other hand, under the null 
charge surface, because of the pair creation, the 
charge density is the combination of the density of 
positrons and electrons, p{r < rnuii) = P- + 
and n_(r) > \p{r)/e\. Wang et al. (2010) show 
that the 77 of the outgoing particles in the outer 
gap of the 7-ray pulsar is around 0.5, therefore, 
we chose 77 — 0.5 for the incoming particles. 

2.3. The outgoing and incoming particles 

The evolution of the Lorentz factor of a particle 
moving along the magnetic field can be described 
as 

where s{r) is the curvature radius of the mag- 
netic field line at r*, the first term on the right 
side is the power of the curvature radiation, and 
the second term is the acceleration of the electric 
field. From the null charge surface to the light 
cylinder, where the accelerating electric field i^y is 
strong, the outgoing particle can reach an equilib- 
rium state such that it radiates the energy it ob- 
tains from the electric field as curvature radiation, 
which is described as eE^\{f)c ~ 2e^C7^(r)/3s^(r), 
by requiring dj/dt = 0. The Lorentz factor of the 
outgoing particle in the equilibrium state is thus 



given by 

and the corresponding characteristic energy of the 
emitted curvature photon is given by 

E .-^ 3 fe7,V(r-) 

For a particle moving from the null charge sur- 
face to the stellar surface, its Lorentz factor is cal- 
culated by solving Equation ([5]) using the Runge- 
Kutta method. Because the electric field in this 
region is much reduced, the energy loss by the cur- 
vature radiation of the particle cannot be compen- 
sated for by the potential drop of the gap. This sit- 
uation has been discussed in Takata et al. (2010). 
When the Lorentz factor of the particle drops low 
enough, the curvature energy loss time scale be- 
comes comparable to the time scale of the parti- 
cle's movement to the stellar surface. There is a 
minimum energy ~ lOOMeV for the curvature pho- 
ton of a particle that is purely losing energy, and 
this energy does not depend on any pulsar parame- 
ters nor the curvature radius of the local magnetic 
field. As will be shown below, the energy of the 
synchrotron radiation is contributed by incoming 
particles between 20i?s and 50i?s. Our numerical 
calculation shows that the characteristic energy of 
the emitted curvature photon Ein (r) in this region 
spans from 200MeV to 400MeV. 

3. The pulse profile of the synchrotron ra- 
diation with a significant pitch angle 

Magnetic pair creation favours larger values of 
6p, the pitch angle between the photon's propaga- 
tion direction and the magnetic field. No magnetic 
pair creation will happen, no matter how strong 
the magnetic field is, if the photon's direction is 
parallel to the magnetic field line, which means 
-B|| = 0. Conservation of momentum requires that 
the secondary pair generated by magnetic pair cre- 
ation makes the same angle to the magnetic field 
as that of the original curvature photon, while the 
opening angle of the pair resulting from pair cre- 
ation can be neglected compared with this pitch 
angle. The condition for magnetic pair creation 
cannot be satisfied at the birth place of the curva- 
ture photon, where the photon moves tangentially 
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to the magnetic field line. The curvature pho- 
ton should move along its path so as to encounter 
enough perpendicular magnetic field, and the sec- 
ondary pair can therefore obtain a significant pitch 
angle from the pair creation, which cannot be ne- 
glected in the calculation of the pulse profile. 

The equations of the motion of the particle in 
the acceleration region are (Harding, Usov & Mus- 
limov 2005; Hirotani 2006; Takata, Chang & Shi- 
bata 2008): 



2e^j^B^ sin^ e„ 



'in 



and 
dP 



dt 2B ds ^ 



26^7252 sin^ I 



■ sm( 



(8) 



(9) 



If the pairs are generated in the acceleration re- 
gion, according to Equation ([5]) , the electric field's 
acceleration can reduce their pitch angles. On the 
other hand, if the pairs are generated in the re- 
gion without strong accelerating electric field, and 
we neglect the first term of the right hand side 
of Equation ([9]), which represents the adiabatic 
change along the dipole field line, the motions of 
the secondary pairs can be described by. 



dP|| 

~dr 



and 



26^72^2 sin^ 



2527252 sin2 ( 



■ cos( 



pi 



(10) 



(11) 



dt irn^d^ 
which indicate that the pitch angle Op stays con- 
stant while the particle loses its energy by emitting 
synchrotron photons. 

In adopting our previous three dimensional 
studies of the radiation of pulsars, we introduce 
a factor a to represent the magnetic field lines 
in the gap and take a = 1 for the last-open 
field lines. We trace the field lines to determine 
the coordinate values [Xo(0p), io(0p), ^o(0p)] 
of the last-open field lines at the stellar sur- 
face. The coordinate values of an arbitrary layer 
of the magnetic field lines with a given "a", 
[Xq(0p), yQ((/)p), Zq((/)p)], can be determined by 
using X'^((j)p) = aXo((?!)p), ^o'l'/'p) = oY^i^'p) and 
Z^((/.p) = {Rl ~ a'X'^ - a2yo'2)i/2_ ^he relation 
between the z and the a is approximated as 

z{x,(j)p) = -^ — —h2{x,(t)p), (12) 



where amin corresponds to the upper boundary of 
the gap. In this paper, the dimensionless thickness 
amin is independent of (j)p and chosen as amin = 
0.925, which corresponds to the fractional size of 
the gap / we choose. 

We divide the region of 0.925 < a < 1 into 30 
layers, which are represented by {i — 1, ...,30). 
We use 180 field lines to divide each layer into 180 
tubes. Because of the inclination of the magnetic 
axis, some tubes' null charge surfaces are too close 
to the light cylinder to form acceleration regions. 
Therefore, in these 180x30 tubes, we consider the 
radiations from the ones with (j)p £ [—70°, 100°], 
where the (pp is the polar angle of the field line on 
the stellar surface. 0p = 0° is at the right-hand 
side of the north pole's magnetic axis in the plane 
including the rotation axis and the magnetic axis, 
where the rotation axis is at the left-hand side of 
the magnetic axis of the north pole. 

We trace the field lines in the gap, from the 
stellar surface to the light cylinder. At each point 
on the field line between the stellar surface and 
null charge surface, we follow these steps: 

1 We calculate the direction of the curvature 
photon Vcur (Takata, Chang & Cheng 2007): 



= XvpB/B + nx r, 



(13) 



where the first term represents the motion 
along the magnetic field line, Vp is calcu- 
lated from the condition that \vcur\ = c, and 
the second term is the drift motion. If the 
magnetic field line is going outwards from 
the stellar surface, the outgoing direction 
is represented by A = 1, and the incom- 
ing direction is represented by A = — 1. If 
the magnetic field line is going towards the 
stellar surface, the outgoing and incoming 
directions are represented by A = — 1 and 
A = 1, respectively. 

2 We trace the Vcur of the curvature photon 
emitted by the incoming particle to find the 
location of the pair creation. The thermal X- 
ray photons of typical temperature of a neu- 
tron star cannot convert the incoming cur- 
vature photons with hundreds of MeV into 
pairs. However, we cannot rule out the pos- 
sibility that the keV soft synchrotron pho- 
tons emitted from the gap of another pole 



7 



can enter the gap and collide with the in- 
coming curvature photons to become pairs, 
therefore we integrate both of the optical 
depths of the magnetic pair creation and 
photon-photon pair creation along the path 
of the curvature photon. 

The optical depth of magnetic pair creation 
for a curvature photon with the energy 
is 

tb{E^)^ Jr\r,E^)dx, (14) 

where l{f,Ej) is the mean free path of the 
curvature photon of E^ at position r, which 
is given by (Erber 1966) 

l{f,E^) = -J^:;^ B±lr) '^^'P^sxi^-E^)')^ 

(15) 

where Bq m?c^/eh = 4.4 x lO^^G, B± = 
B sin 9p, and x ^ 1 is required. Here 
the rotation of the magnetic field is consid- 
ered, which in fact can be neglected for PSR 
B1509, because of its period of 150ms. For 
the photon emitted at i = 0, the magnetic 
field at its position r when t — t' is 

B{(j),t') ^ MB{(j)~m' (16) 

where the (p is the position angle of the pho- 
ton &t t — t' , M is & rotation matrix that 
rotates the magnetic field for Vlt' in the di- 
rection of position angle, and the position of 
the photon at t' and its birth place satisfy 

r - To = t'Vcur- 

The optical depth of photon-photon pair cre- 
ation is 

Tp{f,E^)^ j j nsoft{e,r)a-{E^,e)dedx, 

(17) 

where the cross section for photon-photon 
pair creation between the curvature photon 
with Ej and the soft photon with e, a{E^, e) 
is given by 

a(E„ e)^f-^aT{l-v' m - v') H'^) - - 
V = - {rrieC^Y / E^t, 

(18) 

where ctt is the Thomson cross section. We 
assume that the angle between the direc- 
tions of the two photons is 90° on aver- 
age. The number density of the soft photons. 



"■so/t(e, ^ is estimated as, 

n..,.(.,rO=0.1^J, (19) 

where the Fobs is the observed flux of soft 
photon radiation in units of MeV/cm^/s, 
r is the distance to the center of the star, 
and d is the distance of the pulsar. We use 
the observed data of 2—20 keV, the energy 
range of a photon that can convert pho- 
tons of hundreds of MeV to a pair, provided 
by Ginga (Kawai et al. 1992). The value 
of d is chosen as 4 kpc (Cordes & Lazio 
2002). Not all of the observed soft photons 
can enter the gap of another pole and have 
collisions under large angles with the curva- 
ture photons, therefore, instead of making a 
simulation to trace the paths of these soft 
photons, we assume 10% of them join the 
photon-photon pair creation. However, as 
will be shown in the next section (Figure [5]), 
the pairs emitting synchrotron photons are 
mainly generated by magnetic pair creation, 
and the photon-photon pair creation's con- 
tribution is negligible. 

3 We integrate the two optical depths given 
above along Vcur, and calculate the direc- 
tions of the synchrotron photons, when 
ArB(10MeV< E^ < 1.5GeV) > 0.1 or 
Arp(10MeV< E^ < l.SGeV) > 0.1. These 
conditions mean that 1 - e^'^^^^^-'^ > 0.095 
or 1 — e-^'^p(^"t) > 0.095 of the curvature 
photons of E-y at r become pairs. The Vgyn 
in the hollow cone of the synchrotron radia- 
tion can be expressed as, 

Vsyn = >^VpV',y,J\v',yJ + ^ X f, (20) 

where Vp is calculated from the condition 
\vsy7i\ = c and v'^y„ is given by, 

"syn,x ^ + tan0p{uoo cosT + Vx sinT) 
'^)]) v'gyj^y = By +tan9p{uyC0sT + VySinT) . 
v'gy^ ^ = Bz + tan 9p {uz cos T -I- U2 sin T) 

(21) 

In the equations above, B is the normalized 

local magnetic held, u = ( i^g2^BlY/2 , - (^2+^2)1/2 : 0) 

and V — are two orthogonal vectors 

\Bxu\ " 
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perpendicular to the magnetic field. Op is the 
pitch angle of the pair, which is the same as 
that of the curvature photon. For photons 
of hundreds of McV, the opening angle of 
the pair obtained from the pair creation is 
much smaller than the pitch angle of the 
curvature photon and therefore is negligible. 
If Op < 90°, then A = 1; and if Op > 90°, 
A = — 1. T G [0, 27r) represents a direction 
along the cone, we trace the Vgyn of each T 
and remove the photons that will enter the 
region dominated by the multi-pole field. 
The velocities of the photons will be used to 
calculate their pulse phases, therefore they 
should be measured at the same time, and 
the Vsyn should be rotated — ilt' in the di- 
rection of position angle. 

4 Having the direction of the photon, Vsyn or 
Vcur, the viewing angle Cj which is the an- 
gle to the rotation axis of the emission's di- 
rection, and the pulse phase of the syn- 
chrotron or curvature photon can be ob- 
tained from (Yadigaroglu 1997) 

r cosC = Wz/|w| ^22) 

In the skymap (e.g. Figure [2]), where the 
X-axis is the pulse phase ip and the Y-axis 
is the viewing angle the highly beamed 
curvature photons emitted from r are repre- 
sented by a point C); on the other hand, 
the synchrotron photons emitted from r are 
represented by a line, as the T of Equa- 
tion ([21]) increases from to 27r. In practice, 
the increment of T, the AT, should make 
|AC| and \Aip\ stay nearly constant. A vis- 
ible synchrotron photon satisfies |C — /3| < 
0.5° and its 

"^syn should not point to the re- 
gion dominated by the multi-pole field. 

Figure [2] shows the differences between the 
{ip, C) of curvature radiation and those of the syn- 
chrotron radiation with a significant pitch angle. 
In this figure, the curvature photons (dashed line) 
and the synchrotron photons (solid lines) originate 
from the region between the null charge surface 
and the stellar surface of only one magnetic field 
line. In the direction of the curvature photons 
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Fig. 2. — The skymap of the viewing angle C 
and the pulse phase ip of the synchrotron radia- 
tion (solid line) and the incoming curvature radi- 
ation (dashed line) originating from one magnetic 
field line, where the inclination angle is a = 20°. 
The incoming curvature photons, which are repre- 
sented by the points on the dashed line, are emit- 
ted by the particles moving from the null charge 
surface to the stellar surface of one magnetic field 
line. The synchrotron photons are emitted from 
the place, where ATi3(10MeV< E < 1.5GeV) > 
0.1 or ATp(10MeV< E < 1.5GeV) > 0.1. Along 
the directions of the curvature photons, there are 
many places where these conditions are satisfied. 
The hollow cones of the synchrotron radiation are 
represented by the solid lines. 

emitted from r, which is represented by a point 
on the dashed line, there are about r/Ar hollow 
cones with different pitch angles, where r is the 
optical depth of the pair creation integrated from 
r to the light cylinder along Vcur, and At = 0.1. 
Each hollow cone is represented by a solid line in 
this skymap. Because the hollow cone of the syn- 
chrotron radiation of the secondary pair contains 
the direction of the original curvature photon, the 
lines of the synchrotron radiation in the skymap, 
which have the same birth place of the original 
curvature radiation, intersect at the (ip^ C) of the 
original curvature radiation. A large step length is 
used to trace this field line, hence there are nearly 
five points of the curvature radiations, where the 
lines of synchrotron radiation intersect. 

We use a smaller step length to trace the field 
lines within cpp e [-70°, 100°] of a = 0.925 - 1.0. 
The synchrotron radiation of the secondary pairs 
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is shown by the color map in Figure |3l where the 
inclination angle a = 20°, the black lines and the 
gray lines are the outgoing curvature radiation and 
the incoming curvature radiation emitted from the 
field lines of a = 0.95, respectively. The outgoing 
curvature photons can be missed by the viewing 
angle < 60° under a — 20°. Furthermore, if the 
viewing angle is smaller than the inclination an- 
gle, even the survival incoming curvature photons 
can be missed. Therefore, small inclination angle 
and viewing angle may be the reason for the dis- 
appearance of the GeV photons. In this figure, 
the dot circles are the pulse phases and viewing 
angles of the two polar caps. It is widely accepted 
that the radio emission comes from the polar caps 
for most of the 7-ray spin-down powered pulsars, 
and the pulse phase of the observed radio emis- 
sion is that of one of the two polar caps, which 
is set to be here. Figure |4] lists the light curves 
of two rotation periods of the outgoing curvature 
radiation (dashed lines) and the synchrotron ra- 
diation (solid lines) of the pairs generated by the 
magnetic pair creation. The light curves of the 
two kinds of radiation have different scales in the 
y-axis. As the inclination angle and the viewing 
angle decrease, the pulse profile of the synchrotron 
radiation gets closer to the observed one of PSR 
B1509, i.e. a wide single peak. So, a small inclina- 
tion and a small viewing angle are required both 
by the pulse profile of the soft 7-ray radiation and 
the absence of the GeV photons. 

The light curve of the synchrotron radiation 
read from the color map, however, is close to but 
not the final result, because it is not energy de- 
pendent and the calculation of the color map does 
not consider the local emissivity or the attenua- 
tion of the original incoming curvature photons, 
which can change the shape of the light curve. In 
order to get the energy dependent light curves, we 
have to calculate the phase resolved spectra. 

4. The spectrum and the energy depen- 
dent light curves of the X-ray and soft 
7-ray emissions 

The secondary pairs are generated in the region 
where the magnetic field can satisfy — — ^ 1, 
where exnx is the energy density of the soft pho- 
tons that can collide with the pairs. This means 
that the pairs lose their energy mainly via syn- 
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Fig. 3. — The skymap of the viewing angle ( 
and the pulse phase ijj of the outgoing curva- 
ture photons (black line) emitted from the null 
charge surface to the light cylinder, the incom- 
ing curvature photons (gray line) emitted from the 
null charge surface to the stellar surface, and the 
synchrotron photons (darkness) emitted from the 
place where ArB(10MeV< E < 1.5GeV) > 0.1 or 
Arp(10MeV< E < 1.5GeV) > 0.1. The curvature 
photons are from the lines of a = 0.95, but the 
synchrotron photons' C and ip are calculated by 
tracing the field fines within (/)p € [-70°, 100°] of 
a = 0.925 - 1.0. The inchnation angle is a = 20°. 

chrotron radiation, and Inverse Compton scatter- 
ing can be neglected. Because the observed spec- 
tra are non-thermal, we neglect the blackbody ra- 
diation from the stellar surface. 

We follow the four steps described above to find 
the positions where the radiations' directions sat- 
isfy IC — /3| < 0.5°. The photon spectrum of the 
synchrotron radiation emitted from r is given by 

d^Nsyn{E,f) ^ V3e^Bir)smep /-^"-^ dN,{f) ^ 
dEdt i'nmeC^hE Je ■ dE^ 

(23) 

where the E^ax and £'„„;n are the maximum and 
minimum energy of the secondary positron or elec- 
tron, respectively. The energy of the pair should 
be in the range of mc^ < E^ < i?ctir/2, where 
the Ecur is the energy of the original curvature 
photon. F{k) = K5/3{^)d£_, where is 
the modified Bessel functions of order 5/3, and 
K = E/Esyn{r, Ee), Esy„{f, Ee) is the critical syn- 
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chrotron photon energy, which is given by 



'3heB{r}sin9p{r)E^ 



(24) 



The energy distribution of the secondary pairs at 
r is 



■a. 



dE~, dt 



E,{r,e,) Je 
g-r(i5^,r-)(^ _ g-Ar(£;^.r))^£;/^ 

(25) 

where Ee{r,E^) = 2e'^B'^{r)sm^ ep{r)El/3mlc'' 
is the synchrotron energy-loss rate of the par- 
ticle with energy i?e, Emax is the maximum 
energy of the original curvature photons, and 
T{Ej,r) = TB{Ej,r) + Tp{Ej,r) is the integrated 
optical depth from ro, the place where the orig- 
inal curvature photons are emitted, to r. The 
definitions of the two optical depths are given by 
Equation (HH) and Equation ((HI)- ^T{E^,r) is 
the sum of the two optical depths between r and 
f+ Af, and Ar is the step we take, (r — tq), Ar 
and Vcur are parallel. 

The photon spectrum of the curvature radiation 
emitted from fo, is 



dEjdt 2Trhs{fo)E-y 



F(«), (26) 



where s(ro) is the local curvature radius of the 
magnetic field line, F{k) is the same as that in 
Equation (1^^ . but k = E^ / Ecuriro), while Ecur 
and 7 are the characteristic energy of the ra- 
diated curvature photons and the corresponding 
Lorentz factor of the particle respectively. For 
incoming curvature radiation, E^ur is obtained 
by solving Equation ([5]); while for outgoing cur- 
vature radiation, E^ur is given by Equation ([T]). 
AA'^(ro) = Ay(ro)n(ro), where n(ro) is the num- 
ber density of the accelerated particle, and Ay(ro) 
is the volume of the tube-like region at tq , which 
can be obtained from the conservation of the mag- 
netic flux, 



AV{r^) = 



2T:rph2iRs)B,Al{fo) 
NANgBiro) 



(27) 



where Al is the step length used to trace the field 
line, Tp is the polar radius of the gap, and Na x 
Nb = 180 X 30 is the number of the lines we use 
to divide the region of a = 0.925 — 1. For the 
outgoing particles outside the null charge surface, 



the number density n(ro,ro > Tnuii) = i^5(ro)[l — 
g{z)]/2Trec; for the incoming particles inside the 
null charge surface, the number density is defined 
in Equation 77 is chosen as 77 — 0.5, which 
equals that of the outgoing particles. 

If the incoming curvature radiation is emitted 
far enough from the stellar surface, the survival 
radiation from the attenuation caused by the pair 
creation has the chance to be observed at certain 
viewing angles, given by 



{rQ,E^) _ d'^Ncuri'ro,Ey) 



dE-ydt 



dE~fdt 



exp(-T(i;^,r')), 
(28) 



the r' is the position where we stop tracing the 
path of the curvature radiation and integrating 
the two optical depths along Vcur- Some of the 
synchrotron photons may penetrate to the region 
where the magnetic field can convert them into 
pairs. Therefore, we trace the trajectories of the 
synchrotron radiation with |C — /3| < 0.5° along 
Vsyn and integrate the optical depth of the mag- 
netic pair creation. Here, photon-photon pair cre- 
ation can be neglected. The survival synchrotron 
radiation is given by 



d^N. 



syn,sur 



{r,E) _ d^N,yr.{r,E) 



dEdt 



dEdt 



exjp{-TBiE.f,r')), 
(29) 

where r' is the place where tb becomes constant. 

The total photon flux received at Earth is the 
sum of the survival synchrotron radiation, sur- 
vival incoming curvature radiation and the out- 
going curvature radiation, 



Ftot{E) - -^E 



d^N^yr,,s^r-i,E,ri) 



D^^'i-^i Ansya{ri) dEdt 

1 rf''JVc^r.^,^.s^,^(-E,rj) 



dEdt 



-I- 



dEdt 



(30) 

For PSR B1509, the curvature radiation is outside 
the viewing angle. The distance of PSR B1509 
is chosen as D = 4kpc, which does not conflict 
with the H/ absorption measurement: 5.2±1.4 kpc 
(Gaensler et al. 1999), and the dispersion mea- 
surement: 4.2±0.6 kpc (Cordes & Lazio 2002). 
The local solid angle of the curvature radiation 
is given by 



A17™r(r) ^ nO' 



(31) 



where 9 is chosen as 0.5°, according to the criterion 
of the visible radiation |C — /?! < 0.5°. Because the 
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secondary pairs are generated at places without 
an accelerating electric field, as shown in Equa- 
tion ([T0|) and Equation (fTTj) . their pitch angles are 
constant, and the sohd angle of the synchrotron 
radiation emitted from r is that of a hollow cone, 
which is given by 



AUsyn {r) = 27r 



'P^ 2-1 



(32) 



where 7 is the Lorentz factor of the particle. 
Af2sy„(r) depends on the energy of the particle, 
as the particle loses its energy via synchrotron 
radiation, I/7 becomes larger and the resultant 
AJlsi/n becomes larger. This makes the observed 
flux smaller. However, at the same time, the in- 
visible synchrotron radiation emitted from other 
places becomes visible under larger I/7. We as- 
sume that the two effects cancel out each other, 
and use Arisyn(7 = 250MeV/meC^). 

By fitting the observed soft 7-ray spectrum and 
the energy dependent light curves, we find that 
the inclination angle and the viewing angle of PSR 
B1509 are 20° and 11°, respectively, which are dif- 
ferent from those of the fitting by Romani & Yadi- 
garoglu (1995) and the fitting by Zhang & Cheng 

(2000) , where the inclination angle is chosen as 
a = 60° . The fitting of the observed polarization 
profile of the radio emission, by using the rotating- 
vector-model (Radhakrishnan & Cooke 1969), can 
provide a constraint on the inclination angle. The 
fitting of the PSR B1509 data by Crawford et al. 

(2001) suggests that a < 60°, which does not con- 
flict with our result. 

Figure [5] shows the spectra of a = 20° and 
j3 — 11°, where the outgoing curvature radia- 
tion and the possible survival incoming curvature 
radiation are missed by the line of sight. The 
dashed and solid lines are the cases with and 
without the attenuation of the synchrotron radi- 
ation of the secondary pair caused by the mag- 
netic field, and show this process cannot be ne- 
glected. The reason for the low cut off energy of 
the spectrum of PSR B1509 is that the magnetic 
field converts some of the synchrotron photons 
into pairs. This figure also shows the contribu- 
tion of the pairs generated by photon-photon pair 
creation (dot-dashed line), which confirms that 
the photon-photon pair creation is negligible and 
the pairs emitting synchrotron photons are mainly 
generated by the magnetic pair creation. Because 



these pairs are generated close to the stellar sur- 
face, where the soft photon density is high enough 
to allow photon-photon pair creation, their syn- 
chrotron radiations above IMeV are strongly at- 
tenuated by the magnetic field. 

Figure |6] compares the simulated phase aver- 
aged spectrum with the observed one. The simu- 
lated spectrum has different shapes above and be- 
low 1 MeV. Below IMeV the spectrum is the same 
as that without the attenuation caused by the 
magnetic field; while above IMeV the spectrum 
has an unusual shape with a low cut off energy. 
There is an estimation of the optical depth of the 
magnetic pair creation by Ruderman & Suther- 
land (1975), which is, when x in Equation ([T5|) 
satisfies 

X~' « 15, (33) 

tb is unity. Here we define a parameter Ecri{r), 
which is the maximum energy of the photon that 
can survive from the magnetic pair creation at the 
position f as 



Ecriir) = 



2mc^ Bq 
15 Bi_{r)' 



(34) 



Ecri = IMeV corresponds to the maximum energy 
of synchrotron photons surviving from the mag- 
netic pair creation. The maximum these pho- 
tons encounter is about 3x lO^^G, so the spectrum 
with E < IMeV can maintain its shape of the syn- 
chrotron radiation, but the part with E > IMeV 
is changed by the magnetic field. 

To obtain the energy dependent light curves, we 
bin the photons by their pulse phases and energy. 
The number of the photons with Ei < E < E2 
measured at pulse phases between ijji and '4>2 is 
calculated from 

N^{Ei,E2,^i,^2) (X / Ftot{E, ^1 < < iJ2)dE. 

(35) 

As the secondary particle loses its energy via syn- 
chrotron radiation, its Lorentz factor 7e becomes 
smaller and the hollow cone can be observed for 
longer. In order to take into account this fact, we 
redistribute the photons with energy E oi ip into 
± 7~^, where 7e is the Lorentz factor of the pair 
emitting these photons with typical synchrotron 
energy E under the local B±{f). Figure [7] shows 
the integrated energy dependent light curves in 
the plane of pulse phase and energy, where there 
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Fig. 4. — The light curves of two rotation periods 
of the synchrotron radiation (sohd Unes) ignoring 
the attenuation of the original incoming curvature 
radiation, and the light curves of the outgoing cur- 
vature radiation (dashed lines), of different incli- 
nation angles and viewing angles. The two types 
of radiation have different scales on y-axis. 
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Fig. 5. — The spectra of the synchrotron radia- 
tions of the secondary pairs generated by pair cre- 
ations. The solid line is the one without consid- 
ering the attenuation of the synchrotron photons 
by the magnetic pair creation. The dashed line 
is the spectrum of the survival synchrotron pho- 
tons from magnetic pair creation. The dot-dashed 
line is the contribution of photon-photon pair cre- 
ation, which considers the attenuation caused by 
the magnetic pair creation. 
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Fig. 6. — The simulated phase averaged spectrum 
of PSR B1509, compared with the observed data. 
The triangles are the 2(T upper limits for the to- 
tal fluxes in the EGERT energy domain. The fllled 
circles arc the COMPTEL flux points derived from 
the excess counts in the 0.15—0.65 phase range, 
and the open circle is the 10— 30MeV flux in the 
0.15—0.65 phase range above the spatially deter- 
mined background (Kuiper et al. 1999). The filled 
squares are the observed pulsed fluxes by AGILE 
in the 30-lOOMcV and 100-500MeV bands (Piha 
et al. 2010). The filled diamonds are the 2a upper 
limits of Fermi LAT derived with gtlike, and the 
open diamond is the integrated flux derived from 
the light curves observed by Fermi (Abdo et al. 
2010). The lines are the X-ray spectra observed 
by ASCA (dotted line) (Saito et al. 1997), Ginga 
(dot-dashed line) (Kawaiet al. 1992), OSSE (long- 
dashed line) (Matz et al. 1994), Welcome (dot- 
dot-dashed line) (Gunji ct al. 1994), and RXTE 
(short-dashed line) (Marsden et al. 1997). 
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are 100 bins in the axis of energy and 180 bins in 
the axis of pulse phase, which means that there 
are 100 energy dependent light curves integrated 
from 180 phase resolved spectra. The darkness 
represents the percentage of the number of the 
photons of a certain pulse phase interval, in the 
total number of photons of certain energy range. 
Below lOMeV, there is an asymmetric wide sin- 
gle pulse at -0 ~ 0.25, and above lOMeV, another 
peak shows up at ^ ^ 0. The explanation for such 
a phenomenon will be given in Section 15.21 We 
compare the energy dependent light curves with 
the observed ones in Figure [H 



4 

m 



5. Discussion 

5.1. The lost 7-ray photons 

The magnetosphere of PSR B1509, in fact, can 
also generate GeV 7-ray photons by accelerating 
the charged particles in the outer gap, just as in 
other 7-ray spin-down powered pulsars. However, 
these outgoing curvature photons are missed by 
our line of sight, which happens more easily un- 
der a smaller inclination angle. And because the 
viewing angle is smaller than the small inclina- 
tion angle, even the survival incoming curvature 
photons are totally invisible. This requires us, to- 
gether with the observed wide asymmetric single 
peak, to choose (3 = 11° and a = 20°, which does 
not conflict with the observed polarization of the 
radio emission. 

Figure |9] shows the phase averaged spectra of 
different inclination angles and viewing angles. 
Under a — 20°, as the viewing angle increases, 
the incoming and outgoing curvature radiations 
become visible. In a huge magnetosphere with 
a small inclination angle, the incoming curvature 
photons have a chance to escape from being con- 
verted into pairs by the magnetic field. Along one 
field line, the place of /3 = 70° is closer to the 
null charge surface than that oi /3 = 50°, there- 
fore, the incoming curvature photons have higher 
energy. When a becomes larger, such as 45°, the 
outgoing curvature radiation shows up earlier than 
the incoming one when /3 increases. To obtain a 
clear picture of the connection between the two- 
layer structure and the region emitting incoming 
curvature photons besides the null charge surface, 
a general multiwavelength study of the high en- 
ergy radiations of the pulsars is needed, where the 



Fig. 7. — The simulated energy dependent light 
curves in the pulse phase and energy plane. The 
darkness represents the percentage of the number 
of photons of a certain pulse phase interval, in the 
total number of photons of a certain energy range. 




Fig. 8. — The simulated energy dependent light 
curves (gray histograms), compared with the ob- 
served data (solid lines), which are provided by 
Fermi LAT (Abdo et al. 2010), AGILE (Pilia et 
al. 2010), COMPTEL (Kuiper et al. 1999), OSSE 
(Matz et al. 1994) and BeppoSAX (Cusumano et 
al. 2001). 
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different a and /3 are considered. 

5.2. The second peak of > lOMeV 

The observations (Kuiper et al. 1999; Abdo et 
al. 2010; Pilia et al. 2010) indicate that there 
is another peak around ■0 ~ in the hght curve 
when E > lOMeV. Here the energy dependent 
light curves (Figure [TUl Figure[T31), phase averaged 
spectra (Figure [T2t and the skymaps of the syn- 
chrotron radiation (Figure llip of the regions with 
different distances from the star under a = 20° 
and (3 = 11° are made to explain this second peak. 
The reason that the emission region is separated 
by the distance from the star, is that the dipole 
magnetic field plays an important role in the re- 
sultant shape of the light curve of the synchrotron 
radiation, by determining the pitch angles of the 
secondary pairs, which are generated via the mag- 
netic pair creation. The pitch angle correlates 
with the distance between the place where the pair 
creation happens and the birth place of the origi- 
nal curvature photons. The farther the curvature 
photon goes before being converted to a pair, the 
larger the pitch angle the secondary pair has. Ac- 
cording to Equation p3|) , the distance of the cur- 
vature photon's penetration is mainly determined 
by two factors: the magnetic field along the path 
of the curvature photon, and the energy of the 
photon. For a fixed B± and E^, if the B{r) is 
higher, then the corresponding 9p{r) is smaller, 
which means a shorter distance between the place 
of pair creation and the birth place of the orig- 
inal curvature photon. And under a fixed B(r), 
the magnetic pair creation of the photon with a 
higher energy E-y requires a smaller pitch angle 

To show how the magnetic field B{r) changes 
the shape of the light curve by changing Op{f), 
we make the light curves of four energy bands, 
which are shown in Figure 1101 whose original 
curvature photons are generated in five regions: 
r < 20Rs, 20Rs < r < SOii^, 30Rs < r < 40i?s, 
40 Rs < r < 50 Rs, and r > 50 Rs- This figure 
shows that the X-ray photons' energy originates 
from particles of r < 50Rs, of which the region of 
SORs < r < 50Rs provides the asymmetric wide 
single peak, and the region of 20Rs < r < 30Rs 
provides a background that smears out the little 
peak of ip ^ 0.65. For the soft 7-rays, the two 
components of the light curves oi E > lOMeV: an 



asymmetric wide single peak at ip '-^ 0.25 — 0.65 
and a sharp peak at t/i '-^ 0, correspond to the re- 
gions of 30Rs <r < 50Rs and 20i?^ < r < 30Rs, 
respectively. The magnetic fields with different 
strengths generate distinct pulse profiles of the 
synchrotron radiation of the secondary pairs. 

Why is there a peak at 0? Suppose there 
are some curvature photons emitted by the parti- 
cles moving along the field lines towards the stellar 
surface at the place f, on their path along Vcur 
given by Equation (jl3p . then the photons with 
higher energy become pairs earlier, which means 
the pairs have smaller pitch angles 9p. In Fig- 
ure [21 the cone of the synchrotron radiation with 
a small 9p is a "circle" around the polar cap, while 
the cone of a large 9p is a "line" going across the 
skymap. When (3 < a, such "circles" contribute 
the observed peak at -0 ~ when E > lOMeV, and 
the "lines" make the wide single peak. As shown in 
Figure [TTJ when r < 30Rs, there are such "circles" 
around the polar cap, while when r > 30Rs, there 
are only "lines" , which means these synchrotron 
cones have large pitch angles. This is because the 
magnetic field is stronger at r < 30 Rs, and, as a 
result, the photons with higher energy can become 
pairs under smaller 9p. For r > 30 Rg, without 
enough strength of the magnetic field, most of the 
photons need large pitch angles (< 90°) to satisfy 
the condition of pair creation, so in the skymap, 
there is no component around the polar cap. 

The attenuation of the synchrotron photons 
caused by the magnetic field also plays an impor- 
tant role in making the second peak. For an indi- 
vidual synchrotron cone at r' , some of the emitted 
photons still have enough energy to become pairs 
when they approach the stellar surface. Although 
this process is not determined by the birth place 
of the curvature photons directly, it can be found 
in Figure [TH that the spectra of different r have 
different levels of attenuation of the synchrotron 
radiation. For the case of r < 20 Rg, few of the 
synchrotron photons with E > IMeV survive from 
the magnetic pair creation. As r increases, more 
and more synchrotron photons can survive, when 
r > 40-Rs, the attenuation is already insignificant. 
The spectrum of the survival synchrotron radia- 
tion of a certain range of r has two components: a 
spectrum without experiencing the pair creation, 
and a partially absorbed spectrum. The unab- 
sorbed component cuts off at the maximum en- 
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ergy of the survival photons, which corresponds 
to the maximum magnetic field they experienced 
along their paths Vsyn- This maximum magnetic 
field becomes smaller as the birth place of the cur- 
vature photon r gets further away from the stel- 
lar surface. These spectra are consistent with the 
energy dependent light curves that are shown in 
Figure [131 From these light curves, it is found 
that the synchrotron radiation, whose pulse phase 
is closer to -0 ^ 0, has a higher chance to survive 
from the magnetic pair creation. This is because 
the synchrotron radiation with smaller pitch angle 
is emitted further from the stellar surface, where 
the photons have a higher chance to survive. 

Why does the peak at -0 ~ show up when 
the energy increases? If the particle is made via 
magnetic pair creation, its typical synchrotron en- 
ergy Esyn and the typical energy of the original 
curvature photon E^ur nearly satisfy E^yn oc Ecur, 
which can be obtained by combining Equation (I24p 
and ((33l) . Higher Ecur leads to smaller 9p. There- 
fore, for the curvature photons emitted from f, 
under /3 < a, the observed spectrum of the syn- 
chrotron radiation, whose pulse phase is closer to 
'0 = 0, has a higher cut-off energy. 

To summarize, for a fixed birth place of the cur- 
vature photons r, the photons with higher Ecur 
become pairs earlier with smaller 9p and further 
from the star, which raises the synchrotron pho- 
tons' chances to avoid from becoming pairs. For 
a curvature photon with Ecur, if it is emitted fur- 
ther from the star, the 9p of the secondary pair 
is larger, because of the smaller magnetic field. 
When f3 < a, the pulse phase of the synchrotron 
radiation is closer to when the pitch angle is 
smaller, and closer to 0.5 when the pitch angle is 
larger. 

Therefore, the second peak of PSR B1509 is the 
synchrotron radiation emitted by the pairs with 
smaller pitch angles, which are converted from the 
ones with higher energy of the curvature photons 
emitted from the region of 20i?s < r < 30i?s . The 
viewing angle /? determines the energy where the 
second peak shows up. If /3 < a, the second peak 
shows up at higher energy when /3 gets closer to 
a, because the observed synchrotron radiation is 
emitted by the pairs with smaller pitch angle 9p, 
which were the curvature photons of higher energy. 
The region of r > 30i?s contributes the wide single 
peak of ~ 0.25 — 0.65. Because each synchrotron 



cone sweeps the line of sight twice on average, in 
the light curve, there are two peaks besides ip ~ 
0.5, which overlap each other, making a wide single 
peak. The rotation of the neutron star make this 
single peak asymmetric. 

6. Summary 

The outer gap model predicts that most pairs 
created inside the gap are around the null charge 
surface and the electric field separates the oppo- 
site charges to move in opposite directions. Con- 
sequently, from the null charge surface to the light 
cylinder the outflow radiation is dominant whereas 
from the null charge surface to the star the in- 
flow radiation is dominant. Since the electric field 
decreases rapidly from the null charge surface to 
the star, the incoming radiation flux is of an or- 
der of magnitude weaker than that of the outgoing 
flux. Furthermore most of the incoming high en- 
ergy curvature photons are converted into pairs by 
the strong stellar magnetic field of the canonical 
pulsars. Based on these model features we pro- 
pose a model to calculate the non-thermal X-rays 
and soft 7-rays emitted by secondary pairs pro- 
duced by converting the incoming curvature pho- 
tons emitted by the outer gap's primary charged 
particles. This model is used to explain the ob- 
served X-ray and soft 7-ray spectra and the en- 
ergy dependent light curves of PSR B1509 from 
lOOeV to GeV. We argue that the line of sight 
of PSR B1509 is in the direction of the incoming 
beam instead of the outgoing beam, otherwise a 
characteristic power law with exponential cut-off 
spectrum with a cut-off energy around a few GeV 
should be observed. 

The secondary pairs spiral around the mag- 
netic field and emit synchrotron photons, some of 
these synchrotron photons can also become pairs 
via magnetic pair creation. In order to avoid see- 
ing the outgoing flux and fit the observed multi- 
wavelength light curves we need to choose a = 20° 
and /3 = 11°. The observed spectrum is the 
synchrotron radiation emitted from the pairs pro- 
duced by the magnetic field that converts the ma- 
jor part of the incoming curvature photons. The 
peak position of the observed spectrum naturally 
occurs around one MeV because photons above 
this energy scale can be converted into pairs by 
the strong stellar magnetic field if they get close 
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Fig. 9. — The simulated phase averaged spectra of 
different inchnation angles a and different viewing 
angles /3. The dashed lines are the synchrotron 
radiations of the secondary pairs, the dot-dashed 
lines are the outgoing curvature radiations, the 
dotted lines are the survival incoming curvature 
radiations, and the solid lines are the total spec- 
tra. 



r<20R, 20R,<rS30R, 
^^^^^s " -i ^^M- ^^^^^^^^^ 




0.0 O.Z 0.4 OG Oa 1.0 0.0 O.Z 04 0.6 0.8 1.0 

¥ ¥ 



Fig. 11. — The skymaps of the viewing angle C 
and the pulse phase ip of the synchrotron photons 
emitted from the place where Atb(10McV< E < 
1.5GeV) > 0.1 or Arp(10MeV< E < 1.5GeV) > 
0.1. r is the distance of the place where the origi- 
nal curvature photons are emitted. 




Fig. 10. — The energy dependent light curves of 
different birth places of the original curvature pho- 
tons. 
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Fig. 12. — The spectra of different regions where 
the original curvature photons are emitted. 
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Fig. 13. — The energy dependent light curves of 
different regions where the original curvature pho- 
tons are emitted. The darkness represents the per- 
centage of the number of photons of a certain pulse 
phase interval in the total number of photons of a 
certain energy range. 

enough to the star. 

The magnetic pair creation requires a signifi- 
cant pitch angle of the photon. For the curvature 
photons emitted from a certain place, the mag- 
netic field and their pitch angles increase along 
their path, the ones with higher energy become 
pairs first, which have smaller pitch angles and 
higher cut-off energies of their spectra, while the 
ones with lower energy need to go further inwards 
to obtain large enough perpendicular magnetic 
field, and they become pairs with larger pitch an- 
gles and lower cut-off energies. The photons emit- 
ted by the pairs with larger pitch angles have a 
higher chance to become pairs further on. There- 
fore the differences between the light curves of dif- 
ferent energy bands are due to the different pitch 
angles of the secondary pairs, and the second peak 
appearing a.t E > lOMeV comes from the region 
near the star, where the stronger magnetic field al- 
lows pair creation to happen with a smaller pitch 
angle. 

Although in this paper our model is applied to 
explain the soft 7-ray spectrum and the energy 
dependent light curves from 100 eV to GeV of PSR 
B1509, it can be extended to other pulsars as a 
general framework for studying the non-thermal 



X-ray and soft 7-ray emissions. 
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